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ABSTRACT 

Absorption lines from the molecules 0H + , H 2 + , and H^" have been observed 
in a diffuse molecular cloud along a line of sight near W51 IRS2. We present 
the first chemical analysis that combines the information provided by all three 
of these species. Together, OH + and H 2 + are used to determine the molecular 
hydrogen fraction in the outskirts of the observed cloud, as well as the cosmic-ray 
ionization rate of atomic hydrogen. Hjj~ is used to infer the cosmic-ray ionization 
rate of H 2 in the molecular interior of the cloud, which we find to be (2 = 
(4.8 ± 3.4) x 10~ 16 s _1 . Combining the results from all three species we find 
an efficiency factor — defined as the ratio of the formation rate of OH + to the 
cosmic-ray ionization rate of H — of e = 0.07 ± 0.04, much lower than predicted 
by chemical models. This is an important step in the future use of OH + and 
H 2 + on their own as tracers of the cosmic-ray ionization rate. 
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INTRODUCTION 



In the past decade, Hjj" has widely become regarded as an excellent tracer of the cosmic- 



ray ionization rate in diffuse m olecular clouds. Surveys of in such clouds ( llndriolo et al. 



20071 ; llndriolo fc McCalll 120121 ) have enabled us to find variations in the ionization rate be- 
tween sight lines, and to build up the distribution function of cosmic-ray ionization rates 
in the nearby interstellar medium (ISM). However, observations of Hjj~ are currently lim- 
ited to background sources with L-band magnitudes brighter than about L = 7.5 mag. At 
this cutoff, OB stars are only feasible as background sources to distances of a few kpc, 
meaning that H^ observations are primarily lim i ted to the local spiral arm (observations to 



ward the Galactic center, e.g., iGoto et al.ll2002l 12008 . l201ll ; lOka et al.ll200a iGeballe & Oka 



20101 . use dust-embedded objects). An alternative method for inferring the ionization rate 
utilizes the chemistry associated with t he formation and destruction of OH + and H20 + 
( IGerin et al.l 1201 Obi ; iNeufeld et al.l 120101 ) , thought to be dependent primarily on hydrogen 
abstractio n reactions with H? and dissociative reco mbination with electr ons. The HIFI in- 
strument (Ide Graauw et al.l 120101 ) aboard Herschel (IPilbratt et al.ll2010l ) has provided the 
first opportunity to observe both OH + and H20 + with very high spectral resolution, thus 
allowing the use of these ions in constraining the cosmic-ray ionization rate. Background 
sources bright enough for THz spectroscopy are widely distributed throughout the Galaxy, 
and targets from the PRISMAS^ key programme range in distance from about 1 kpc to 
12 kpc. However, the ionization rate inferred from the oxygen chemistry is dependent upon 
an efficiency factor, e, at which atomic hydrogen ionized by cosmic rays will eventually be 
converted into OH + . In order to determine e, we present observations of OH + , H20 + , and 
Hjj~ in sight lines toward W51 and compare the ionization rates inferred separately from the 
hydrogen chemistry and oxygen chemistry. 



1.1. Hydrogen Chemistry 

The interstellar chemistry of Hjj" is rather simple, and the reactions surrounding this 
molecule are given in the top portion of Table [TJ H^ is formed in a two-step process, 
beginning with the ionization of H 2 by cosmic rays, and quickly followed by a reaction of H^ 
with H 2 . Some H^ is destroyed by dissociative recombination with electrons or by charge 
transfer to atomic hydrogen, but these reactions are generally slow compared to the H^ + H 2 
process^] Cosmic-ray ionization is the rate-limiting step in this process as it is many orders of 



7 PRobing InterStellar Molecules with Absorption line Studies 

8 This is no longer true for the + H reaction at low molecular fraction. 
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magnitude slower than proton transfer from FI 2 to H2, and can be taken as the formation rate 
of H3. The primary destruction mechanisms for Hjj" are dependent on the environment under 
consideration. In diffuse molecular clouds, H3" is predominantly destroyed via dissociative 
recombination with electrons. In dense clouds, however, where the electron fraction is much 
lower, Hjj" is destroyed by proton transfer to neutrals such as CO and O. 



1.2. Oxygen Chemistry 

Reactions involved in the chemistry surrounding OH + and H 2 + are presented in the 
bottom portion of Table [TJ The formation of OH + begins with the ionization of atomic 
hydrogen by cosmic rays. This is followed by endothermic charge transfer to oxygen to form 
+ — a process hig hly dependent upon the relative populations in the fine structure levels 



of atomic oxygen ( jStancil et al.l Il999f ) — and hydrogen abstraction from H2 to form OH + . 
OH + is either destroyed by further hydrogen abstraction to form H 2 + , or by dissociative 
recombination with electrons. The same is true for H 2 + , but H 3 + is only destroyed 
by dissociative recombination with electrons. A steady-state analysis of these reactions is 
employed in Section H] in inferring the ionization rate of atomic hydrogen and molecular 
hydrogen fraction from OH + and H 2 + abundances. 

An alternative means of forming OH + is the reaction of O with ILt. This process 
requires a high molecular hydrogen fraction — such that Hjj~ is formed efficiently from cosmic- 
ray ionization of H2 — and a low electron fraction — such that H^" is predominantly destroyed 
by proton transfer to O, forming OH + . As we will show in Section H] that the OH + and 
H 2 + probed by our observations reside in gas with a low molecular hydrogen fraction, we 
omit this pathway from our analysis. 



2. METHODS 

2.1. Target Characteristics 

Observations in the near to mid infrared at UKIRT (United Kingdom Infrared Tele- 
scope) and Gemini South were made using the embedded cluster W51 IRS2 (a=19 /l 23 m 40.0 s , 
5=+14°31'06"; J2000.0) as a background source. THz observations of W51 in the PRISMAS 



key program were pointed toward the 3.6 cm continuum peaks ei and e 2 ( iMehringerl Il994j ) 
at ot=19 h 23 TO 43.9*, 5=+14°30'30 / .'5 (J2000.0), meaning there is a 65" separation between the 
IR and THz pointings. 
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Previous studies of the W51 star forming complex and giant molecular cloud (e.g., 



Carpenter &: Sanderslll998l ; l0kumura et al.ll200ll ; iBieging et al.ll2010l ; lKang et al.ll201Q[ ) show 



emission in CO at velocities between about 50 km s -1 and 70 km s~ l . This gas is thought 
to be dense and associated with the W51 region. However, ICarpenter fc Sanders! (119981 ) 
also report a cloud in the solar neighborhood at 7 km s _1 , and some weaker features at 15- 
25 km s _1 , which are thought to arise from a diffuse molecular cloud and a primarily atomic 
cloud, respective^ (jNeufeld et al. 2002 ; Sonnentrucker et al. 2010 ). All of these velocity 
components are also seen in H I 21 cm observations (IKod 119971 . see G49.5-0.4e spectrum). 
The 25 km s _1 component is much stronger (with respect to other components) in H than 
CO, supporting the conjecture that this is a primarily atomic cloud. In the present study 
we are concerned primarily with the conditions of diffuse molecular clouds, and so will focus 
mainly on the absorption at 5-7 km s~ l . This material shows a few closely spaced absorption 
components (jSonnentrucker et al.l l2010l ; iGodard et al.l 120101 ) and is est imated to be at a 
distance of about 500 pc using a simple Galactic rotation curve analysis (IGerin et al.l 12011 ). 



but m ay be as close as 100-200 pc based on maps of the nearby neutral ISM ( iWelsh et al. 



2010|). At 500 pc the on-sky separation of 65" between the IR and THz pointings corresponds 



to a physical separation of 0.16 pc, so all observations of the diffuse molecular cloud(s) of 
interest should probe roughly the same material. 



2.2. Observations 



Observations made at UKIRT utilized CGS4 (Cooled Grating Spectrometer 4; iMountain et al 

19901 ) with its echelle grating, ~0f!4 wide slit, long camera, and 3x2 pixel sampling mode 
to yield a resolving power of about 37,000 (resolution ~8 km s _1 ), in combination with a 
circular variable filter (CVF) to select the correct order. Targets were nodded along the slit 
in an ABBA pattern to facilitate the removal of sky background. The R(l, 1)" and R(l, 0) 
transitions of H3 at 3.668083 /im and 3.668516 /zm, respectively, were targeted toward W51 
IRS2 and a Lyr (observed for the purpose of removing telluric lines from the science tar- 
get spectrum) on 2001 May 26. Total integration time on the science target was 33.6 min. 
Observations targeting the -R(O) through R(3) transitions of the v = 1-0 band of 12 CO near 
4.64 /im toward W51 IRS2 and the telluric standard a Aql were made on 2001 May 28 with 
an integration time of 9.6 min on the science target. 



At Gemini South the Phoenix spectrometer (IHinkle et al.ll2003l ) was used with its echelle 



9 These clouds are not reported by Okumura et al. ( 2001 ). Bieging et al. ( 201ol) . and Kang et al. ( 2010h 
because they are outside of the covered velocity ranges. 
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grating and 0'.'17 slit to produce a resolving power of about 70,000 (resolution ~5 km s _1 ). 
The L2734 filter was employed to select the order containing the R(l, I) 1 transition of H^ 
at 3.715479 fim. Observations targeting this transition toward W51 IRS2 and the telluric 
standard a Aql were made on 2010 Jul 23 with a total integration time of 4 min on the 
science target. 

The HIFI spectrometer aboard Herschel was used to observe the N = 1 — 0, J = 
2 — 1, transition of OH + and the N KaKc = l n — o, J = 3/2 — 1/2 transition of ortho - 



H 2 + , for which the strongest hyperfin e components are at 971.804 GHz (iMiiller et al.ll2005[ ) 



and 1115.204 GHz ( IMiirtz et al.l Il998f ) respectively. The observations were carried out in 
dual beam switch (DBS) mode, with the reference beams located 3' on either side of the 
source. We used multiple local oscillator (LO) frequencies, separated by a small offset, 
to confirm the assignment of any observed spectral feature to either the upper or lower 
sideband of the (double side band) HIFI receivers. For the H20 + line, the observations were 
performed on 2010 Oct 29 with 3 separate LO settings in the lower sideband of mixer band 
5a (Observation Identifications [ObsIDs] 1342207693, 1342207694, and 1342207695). For the 
OH + line, observations were carried out with three separate LO settings in the lower sideband 
of mixer band 4a on 2010 Oct 28 (ObsIDs 1342207642, 1342207643, and 1342207644). 



2.3. Data Reduction 



Starting from raw data frames, our reduction procedure for IR data utilized standarad 
IRAF0 routines commonly used in spectroscopic data reduction. Upon extracting one- 
dimensional spectra, da ta were trans fered to Igor where we have macros written to 
complete the reduction ( ]McCallll200ll ). A full descrip tion of the data reduction procedure- 
applicable to both Hjj" and CO data — is presented in llndriolol (120111 ) . 



Herschel data were processed to Level 2 in HIPEo using the standard HIFI pipeline, 
providing fully calibrated spectra with the intensities expressed as antenna temperature and 
the frequencies in the frame of the Local Standard of Rest (LSR). For each of the target 
lines, the signals measured in the two orthogonal polarizations were in excellent agreement, 
as were spectra obtained at the various LO settings when assigned to the expected sideband. 
We combined the data from the multiple LO settings, and from both polarizations, to obtain 



10 http://iraf.noao.edu/ 

11 http://www. wavemetrics.com/ 

12 Herschel Interactive Processing Environment 
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an average spectrum for each line. 



RESULTS 



Reduced spectra are shown in Figure [U All species show multiple absorption features 
resulting from multiple gas clouds along the line of sight. Hjj" and CO have components 
at about 5 km s -1 , 50 km s _1 , and 65 km s _1 , and the 5 km s _1 component corresponds 
to the diffuse cloud of interest in our study. OH + and H 2 + show absorption over a wider 
range of velocities, but the analysis of these spectra is complicated due to hyperfme splitting. 
The green curves show the absorption due to only the strongest hyperfme component (found 
as discussed in Section l-O]) . revealing that OH + and H 2 + exhibit absorption at the same 
velocities as Hj~ and CO. Additionally, OH + shows absorption at 25 km s _1 where HF has 
been observed (ISonnentrucker et al.l 120101). While ou r fit t o the H 2 + spectrum does not 



require a component at 25 km s -1 , IWyrowski et al.l (120101 ) reported H 2 + absoprtion at 
22.5 km s _1 toward W51 at a position 60" away from the PRISMAS pointing (17" from 
the IR pointing). Most likely this difference is caused by the the lower signal-to- noise ratio 
in the PRISMAS spectrum compared to the WISH (Water In Star-forming regions with 
Herschet) spectrum. Note, however, that for OH + and H 2 + the many components used 
to fit absorption features do not necessarily correspond to physical clouds, but are simply 
used as a means to quantify the optical depth as a function of velocity. A more complete 
description of the spectra in Figure [1] is given in the figure caption. 



4. ANALYSIS 
4.1. Extraction of Column Densities 

Absorpti on features due to Hj~ w ere fit with Gaussian functions using the procedure 



described in llndriolo fc McCalll ( 120121 ) in order to determine equivalent widths, velocity 
FWHM, and interstellar gas velocities. All three of the absorption features in the R(l, l) 1 
spectrum were fit simultaneously, and the resulting FWHM and gas velocities were used to 
aid in the simultaneous fitting of the six absorption features (three R(l, l) u and three R(l, 0) 
lines) in the R(l, l) u and R(l, 0) spectrum. Absorption line parameters and column densities 
determined from this analysis are reported in Table [2J Note that the severe blending of the 
R(l,l) u line at ~ 50 km s _1 with the R(1,0) line at ~ 5 km s _1 introduces additional 
uncertainty to the parameters extracted for these lines. 



In the present study we are only concerned with the diffuse foreground cloud, and the 
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feature of interest in the CO spectra is the relatively weak absorption near 5 km s _1 that 
is seen in the R(0) and R(l) lines. In the bottom panel of Figure [1] — noticeable in the 
R(0), R(2), and R(3) spectra — it is seen that the continuum level slopes downward from 
about —10 km s _1 to 35 km s _1 . It is unclear whether this feature is due to a broad 
outflow component or some artifact of the instrument or data reduction, but it hinders the 
determination of the continuum level across the absorption feature of interest, and thus the 
determination of the equivalent width. To remove the contribution to absorption from the 
dense line-of-sight gas, remove the sloping continuum level, and extract the equivalent width 
from the feature of interest, we fit each spectrum with the sum of five Gaussian components. 
Three narrow components between about 30 km s _1 and 80 km s _1 fit the dense cloud 
absorption; one broad component centered near 20 km s" 1 with a FWHM of ~ 60 km s" 1 
removes the sloping continuum level; and one narrow component fits the absorption from the 
diffuse cloud at 5 km s _1 . Extracted line parameters and column densities in the optically 
thin approximation for the diffuse cloud are given in Table |2j In the case of the R(2) and R(3) 
lines only the first four components are used in the fit, and upper limits on the equivalent 
width are determined from the standard deviation on the continuum level after dividing the 
spectrum by the fit. 

The spectra of OH + and H20 + were fit using a procedure similar to that described 



by iNeufeld et al.l (120101 ) . Multiple Gaussian components with adjustable line-center optical 
depths, velocity dispersions, and line centroids were convolved with the hyperfine structure 
for each transition to obtain an optimal fit to the observed spectra. From this analysis, we 
find dN/dv (column density per unit velocity interval) as a function of LSR velocity, such 
that integrating the function between any two velocities provides the column density in that 
interval. Results from this procedure are presented in Table [3j 

Total column densities (i.e., the sum of column densities over all rotational levels) in 
the diffuse foreground cloud at ~ 5 km s -1 for molecules studied in this paper, along with 
species reported in the literature that are necessary for our analysis, are given in Table HI 



4.2. Cloud Conditions from CO 

Observations of CO are useful in constraining physical conditions within a cloud. The 
ratio of column densities in the J = 1 and J = levels gives an excitation temperature of 
4.3 K, much lower than the expected kinetic temperature of ~ 70 K, but higher than the 
temperature of the cosmic microwave background radiation. This excitation temperature 
suggests that both radiative and collisional excitation play a role in exciting the J = 1 state. 
The observed excitation temperature for J = and J = 1 implies that the density of the 
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collision partners H2 and H must be significantly below the critical density for the 7=1—7-0 
transition of CO. An analysis of the level populations that i ncludes radiative and collisional 
excitation along with optical depth effects (based on that in lNeufeld et al.lll995[ ) suggests a 
density of n H ~ 1 00 cm" 3 (where w H = n(H) + 2n(H 2 )), in good agreement with the range 
of values found by iGodard et al.l ( 120101 ) . 

The fraction of carbon in the form of CO is also important in understanding the chemical 
conditions in the ISM. Table H] shows that the relative abun dance of CO with respect to total 
hydrogen is 8.01 x 10~ 7 , while that of C + is 1.14 x 10~ 4 (iGerinl |2012| . in press), implying 
that carbon is primarily in ionized form. If electrons are predominantly the result of singly 
ionized carbon, then the fractional electron abundance, x e = n e /ri}i, can be approximated 
by a;(C + ). Because C + ions (and thus electrons) are 140 times more abundant than CO, 
and because the rate coefficient for proton transfer from H^" to CO is about 100 times slower 
than that for dissociative recombination of with electrons, destruction of Hjj" by CO is 
negligible. 



4.3. Ionization Rate from Ho 



The standard steady-state analysis for the formation and destruction of H3 in diffuse 
clouds gives 



C 2 n(H 2 ) = A;(H3 + |e-HH^ 



T. 



where (2 is the ionization rate of H2, n(X) is the number density of species X, n e is the 
electron density, and k(X\Y) is the rate coefficient for the reaction between species X and 
Y. Equation ([T]) can be re-arranged and substitutions made such that 



C2 = k{Rj\e )x e n }i 



mi 

7V(H 2 



(2) 



as described in llndriolo fc McCalll ( 120121 ). where iV(X) is the column density of species X. 
Variables on the right-hand side of equation (j2j) are determined as follows. The molecular 



hydrogen column density, iV(H 2 ), is estimated from 7V(CH) (jGerin et al.ll2010a, see Table 3] 



herein ) using the empirical relationship A^(CH)/iV(H 2 ) = 3. 5j^; 4 x 10 8 from 



Sheffer et al.l 



(120081 ). The H^-electron recombinatio n rate coefficient, fc(Hj~|e ), has been measu red in 



McCall et al. 2004; Kreckel et al. 



2005 



McCall et al 



20101) with 



(120041 ) shown 



multiple laboratory experiments (e.g. 

consistent results, and we adopt the analytical expression from 
in Table [TJ As mentioned above, a hydrogen density of nn = 100 cm -3 is adopted from 
the CO analysis, and the electron fraction is approximated by x(C + ) = 1.14 x 10~ 4 . Note, 
however, that this may underestimate x e in regions where the cosmic-ray ionization rate 
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is high enough that the abundance of H + rivals that of C + - An estimate of where this 
occurs is given by the model chemistry in iHollenbach et al.l (120121 ) . In their model with 
wh = 100 cm -3 and (2 — 4.6 x 10~ 16 s _1 the resulting electron fraction is about 2x(C + ) in 
regions of low molecular fraction (/h 2 < 0.2), but comparable to x(C + ) in regions of higher 
molecular fraction where Hjj~ is expected to primarily form. As such, we omit this effect 
from our analysis of the Hjj" chemistry. Lastly, Af(tht) is determined from the observations 
presented herein, and is given in Table HI Using all of thes e values and assuming T = 70 K — 
an average spin temperature found from H 2 observations (ISavage et al.lll977t iRachford et al. 
2OO2J, I2OO9J) — we find a cosmic-ray ionization rate of (2 = (4.8 ±3.4) x 10 



1 6 s 1 in the diffuse 



cloud component at 5 km s 1 . The uncertainty in this value is primarily due to the scatter 
in the CH/H2 relation, and an assumed 50% uncertainty in the density. 



4.4. Molecular Hydrogen Fraction from OH + and H 2 + 

A steady state analysis of the H 2 + abundance gives the equation 

n(OH+)n(H 2 )£;(OH+|H 2 ) = n(H 2 0+)[n(H 2 )A;(H 2 + |H 2 ) + n e A;(H 2 + |e-)]. 



(3) 



This ca n be re-arranged to prod uce the abundance ratio relation presented in iGerin et al. 
feplObh and lNeufeld et al.l feoioh : 



n(OH + ) 

n(H 2 + ) ~ A;(OH+|H 2 ) ' fa A;(OH + |H 2 ) 



A;(H 2 + |H 2 ) 2x e fc(H 2 + |e" 



(4) 



where fa = 2n(H 2 )/n H . From this equation, the molecular hydrogen fraction, fa, can be 
determined using the abundances of OH + , H 2 + , and electrons, and the relevant reaction 
rate coefficients as 



H 2 



2x e A;(H 2 + 


e')/k(OB. + 


H 2 ) 


A^(OH+)/A^(H 2 0+) - 


- k(R 2 + 


H 2 )/A;(OH + 


H 2 ) 



(5) 



assuming constant densities and temperature in the region probed. Given the OH + and 
H 2 + column densities found from our HIFI observations and the rele yant rate coefficient s 
from Table [U assuming T = 100 K (H I spin temperature adopted by lGodard et al.ll2010l ). 
we find a molecular hydrogen fraction of /h 2 = 0.04 ± 0.01 . T his is compa rable to the 



low m olecular fractions inferred toward W49N (INeufeld et al. 



20101 ) and W31C (IGerin et al. 



2010b! ) from similar observations. However, this is much lower than the cloud- averaged 



13 If x e — 2x(C + ) in atomic regions as suggested by the 



Hollenbach et al. 



mm models, / Ha = 0.08. 
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molecular fraction of fg 2 = 2iV(H 2 )/(iV(H) + 2iV(H 2 )) = 0.60 found using the H and H 2 
column densities reported in Table HI As such, OH + and H 2 + must reside in the primarily 
atomic, outer layers of the observed cloud(s). 



4.5. Ionization Rate from OH + and H 2 + 

Steady-state chemistry for OH + is given by 

e( H n(R) = n(OH + )[n(H 2 )A;(OH + |H 2 ) +n e fc(OH + |e")], 



(6) 



where Ch is the cosmic-ray ionization rate of atomic hydrogen. Here, the destruction of 
OH + is thought to include all important reactions, but not every H + formed by cosmic-ray 
ionization will even tually lead to OH + . To accommodate this fact, we introduce an efficiency 
factor, e, following iNeufeld et al.l (120101 ). The quantity e is defined as the ratio of the OH + 
formation rate to the cosmic-ray ionization rate of H. Solving for the product e(u yields the 
equation 

iV(OH- 



N(R) 



-n H 



|H 2 ) +x e A;(OH + |e- 



(7) 



Because the OH + is thought to reside in the outskirts of the cloud, while the Hjj~ resides 
in the molecular interior, we use a simple pressure balance argument to estimate % in the 
atomic gas. For a purely atomic exterior and a purely molecular interior, pressure balance 
requires the relation T a n a = T m n m , where the number density of collision partners in the 

n H ,a), and 



n„ 



atomic gas is equal to the hydrogen nucleon density in the atomic region 
the number density of collision partners in the molecular gas is equal to one-half of the 
hydrogen nucleon density in the molecular region {n m = nu >m /2). Taking T a = 100 K, 
T m = 70 K, and riH,m = 100 cm -3 , we find n^ a = 35 cm -3 . This value may be used 
for ?2h in equation ([7]), along with the relevant column densities and rate coefficients, to 
obtain eCtr = (0.21 ± 0.1 1) x 10~ 16 s _1 . Using the scaling between ( 2 and Ch given by 
Glassgold fc Langerl (119741 ). 1.5£ 2 = 2.3Ch, and taking the value of ( 2 inferred from Hg and 
e^H inferred from OH + and H 2 + , we find an efficiency factor of e = 0.07 ± 0.04 for the 
production of OH + via cosmic-ray ionization of atomic hydrogen. By combining equations 
(F2]), ©, and (|7]) the efficiency factor can be written as 



2.3 T m iV(OH + ) iV(H 2 



1.5 2T a JV(H) A^(H+) A;(H+|e- 



■ x 



fc(H 2 0+ 


o 


N(OH+)/iV(H 2 0+) - A;(H 2 + 


H 2 )/A;(OH + 


H 2 ) 



+ A:(OH+|e" 
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demonstrating that e is independent of x e and tin. Note, h owever, that if x e differs between 
diffuse atomic and diffuse molecular gas as suggested by iHollenbach et al.l (120121 ) , then a 
scaling factor corresponding to x eia /x e;m must be added to equation (jSj). As mentioned above, 
this effect may result in x eta /x e ^ m ~ 2, in which case e would increase to 0.14. Given the 
dependence of the dissociative recombination rate coefficients of rhj~, OH + , and H2O 4 " 
the efficiency factor scales as e oc {T a /T m )~ lb for different ratios between the temperatures 
in the atomic and molecular gas. 



5. DISCUSSION 



The efficiency factor we determi ne is much lower than that predicted by PDR mod - 



els computed using the Meudon code (ILe Petit et al. 



where 0.5 < e < 1.0 (see discussion in Neufeld et al. 



2006 



Goicoechea fc Le Bourlotl 120071 ) . 



20101 ). In those models, the chain of 



reactions leading from H + to OH + is broken by recombination of H + or + with electrons, 
both of which decrease e. Although recombination of + is not important, recombination of 
H + , while slo\\0, can compete with the + + H2 reaction at low molecular fraction where 
the + + H — > O + H + back- reaction dominates the reaction with H2 that forms OH + . As a 
result, the reaction network cycles between H + and + , sometimes forming OH + and some- 
times forming H. Because of this mechanism, at /h 2 = 0.04 the efficiency factor is about 
0.5. 

Another property of the Meudon models that can reduce e is the inclusion of state- 
specific rate coefficients for the H + + O reaction, and the relative populations in the fine- 
structure levels of ( 3 Pj)0. The ground (J = 2) level is likely the most populated state, and 
the H + + ( 3 P 2 )0 reaction is dra matically slower at l ow (< 100 K) temperatures than reac- 
tions involving ( 3 Pi)0 or ( 3 Po)0 (jStancil et al.lll999l ). If charge transfer to + is inhibited, 
then H + has more time to recombine with electrons, thus decreasing the efficiency at which 
OH + forms. 

Still, recombination of H + alone cannot account for the small value of e = 0.07 that 
we find. In order to reach lower values of e under diffuse cloud conditions, something other 
than electrons must be removing H + from the gas phase. One possible m echanism for doing 
just thi s is the neutralization of H + on small grains or PAH's (lLisztll2003l ). Recent modeling 
efforts (IHollenbach et al.l 120121 ) that account for grain/PAH neutralization find values of 
e ~ 0.1-0.3, much closer to our observationally-derived value. We can apply these low 
efficiency factors to update ionization rates inferred in previous studies of OH + and H^O" 1 ". 



fc(H + |e - ) = 3.5 x l(T 12 (T/300)- a75 cm 3 s" 1 ; UDFA06 
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Rescaling the lower limit on the ionization rate reported in iGerin et al.l ( I2010bl ) to account 
for e = 0.07, we find Ch > 2.6 x 10~ 18 n(H) s" 1 for th e line of sight t oward W31C. Doing 
the same for the range of ionization rates reported in iNeufeld et al.l (120101 ) toward W49N 
results in 8.6 x 10~ 16 s" 1 < Ch < 17 x 10 



" 16 s" 1 . 



These values are high compared to the 
distribution of ionizatio n rates found using Ht (Ch = (2.3t|(J) x 10~ 16 s _1 , converted from 



the mean value of £2 in llndriolo fc McCallll2012l ). However, this scaling procedure is highly 
uncertain at present, and must be improved by determining e in more cloud components. 

It is also possible that our analysis underestimates e given uncertainty in the rate coeffi- 
cient for dissociative recombination of OH + with electrons. As shown in Table (H this process 
is 10 times slower than all of the other dissociative recombination reactions. However, it has 
not been measured at temperatures relevant to diffuse interstellar clouds using vibrationally 
cold OH + molecules. Under such conditions, it is possible that resonance structure in the 
low-energy collision cross section may increase the l ow-temperature rat e of this reaction. 
Indeed, some resonance structure has been observed (lAmitay et al.lll996l ). but that experi- 
ment did not determine the cross section on an absolute scale. Future measurements of the 
OH + + e~ dissociative recombination cross section using storage ring facilities are urgently 
needed0 

Additionally, it would be advantageous to employ a chemical model that is more com- 
plete than the analytical expressions used herein. Several species have been observed in 
the diffuse molecular cloud toward W51, all of which can be used in constr aining the ambi- 



ent p hysical conditions. However, some species — e.g., H (lKoolll997l ). CH + ( Falgarone et al. 
201 0]), OH + , H2O 4 " — are thought to reside in prim arily a tomic gas, while others — e.g., CH 



(IGerin et al.ll2010ah . HF dSonnentmcker et al.ll2010h . H 2 (INeufeld et al.ll2002l : ISonnentrucker et al 
20101 ) . HCO + , HCN ( iGodard et al.ll2010l ). CO — prefer molecular environments. It seems ap- 
parent then, that even with similar velocity profiles, not all of the observed species are 
spatially co-located. Instead, these atoms and molecules are likely probing different portions 
of a cloud complex, including diffuse atomic outer layers, and a diffuse molecular interior. 
Any chemical model attempting to reproduce the observed abundances along this sight line 
must account for these effects. 

Our analysis of OH + and H 2 + shows that both species must reside in gas of low 
molecular hydrogen fraction (/h 2 — 0.04) in order to explain the observed abundance ratio 
between the two species. Hjj~, however, is not efficiently formed in gas with low molecular 



15 If /c(OH + |e ) is 10 times larger than the value adopted in this study (see Table [lj, then the value of e 
required to bring the cosmic-ray ionization rates inferred from H3" and OH + into agreement increases from 
0.07 to 0.23 
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fractions where the H 2 + H reaction competes with the H 2 + H2 reaction, meaning that 
H^j~ must primarily reside in regions of higher molecular hydrogen fraction. As such, it seems 
likely that there is little overlap in the gas probed by OH + and H 2 + , and that probed by 
H;j~ in diffuse molecular clouds (this changes for gas with a low electron fraction where the 
Hjj~ + O reaction becomes important). 



6. SUMMARY 

We have made observations of Hjj~, OH + , H 2 + , and CO in a diffuse molecular cloud 
along closely spaced sight lines toward W51 IRS2. The cosmic-ray ionization rate of H 2 
inferred from the Hjj~ column density is (2 = (4.8 ± 3.4) x 1CT 16 s _1 . Observed OH + and 
H 2 + abundances yield an estimated molecular hydrogen fraction of /h 2 = 0.04 ± 0.01 and 
a product of e£n = ( 0.21 ± 0.11) x 10~ 16 s -1 in the atomic outskirts of the cloud, where 



e is defined following iNeufeld et al.l (120101 ) as the ratio of the OH + production rate to the 
cosmic-ray ionization rate of H. Combining both ionization rates, we find e = 0.07 ± 0.04, 
such that only 7% of H + formed by cosmic-ray ionization goes on to eventually form OH + . A 
possible explanation for the low O H + formation effic i ency is the neutralization of H + on small 



grains and PAHs, as suggested by lHollenbach et al.l (120121 ). Detailed chemical modeling that 



accounts for this process, the change from atomic to molecular gas with cloud depth, and 
all of the species observed in this diffuse molecular cloud, should provide further insight 
regarding the extent to which H^" coexists with OH + and H 2 + . 
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Fig. 1. — Spectra toward W51 in velocity space. First panel: H^" spectra showing the 
R(l, l) 1 line (top spectrum) and the R(l, l) u and R(1,0) lines (bottom spectrum). Vertical 
lines mark the absorption features at about 5 km s _1 , 50 km s _1 , and 65 km s _1 . In the 
bottom spectrum the R(l, 0) line is set to zero velocity, such that the R(l, l) u features are 
shifted —35.4 km s _1 from LSR velocity. Second panel: Spectrum showing the Nx a K c — 
In — Ooo, J = 3/2 — 1/2 transition of ortho-H20 + . A stick diagram centered at 6 km s _1 
shows the hyperfme structure of the observed transition. The black curve is the observed 
spectrum, the red curve is the fit to that spectrum, and the green curve shows only the 
strongest hyperfme component portion of the fit. Third panel: Spectrum showing the 
N = 1 — 0, J = 2 — 1 transition of OH + . Colors are the same as for H20 + . Fourth panel: 
Spectra of the v—1-0 band of 12 CO showing (from top to bottom) the R(3), R{2), -R(l), and 
R(0) transitions. 
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Table 1:: Reaction Network 







Rate Coefficient 




Reaction 




(cm 3 s _1 ) 


Reference 


H 2 + CR - 


> H+ + e~ + CR' 


C2 




H+ + H 2 - 


>H+ + H 


A;(H+|H 2 ) = 2.08 x 1(T 9 


1 


H+ + e-- 


>H 2 + HorH + H + H 


k{H£\e-) = -1.3 x 10" 8 + 1.27 x lO^T" - 48 


2 


H + CR^ 


H+ + e~ + CR' 


Ch 




H+ + -> 


0+ + H 


jfe(H+|0) = 7.31 x KT 10 (T/300) - 23 exp(-225.9/T)t 


3,4 


0+ + H 2 - 


> 0H+ + H 


fc(0+|H 2 ) = 1.7 x 10- 9 


5 


0H+ + H 2 


-> H 2 0+ + H 


fc(OH+|H 2 ) = 1.01 x 10~ 9 


6 


0H+ + e" 


— >■ products 


&(OH + e~) = 3.75 x 10- 8 (T/300)-°- 5 


7 


H 2 0+ + H 


2 ^H 3 0+ + H 


fe(H 2 0+|H 2 ) = 6.4 x 10~ 10 


8 


H 2 + + e" 


— > products 


fe(H 2 + |e-) = 4.3 x 10- 7 (T/300)-°- 5 


9 


H 3 0+ + e" 


— > products 


A;(H 3 + |e-) = 4.3 x 10- 7 (T/300)-°' 5 


10 



Notes: Both £ 2 and Cu have units of s . Rate coefficients are those used in UDFA06 
(www.udfa.net). All temperature dependent coefficients are in terms of the gas kinetic 
temperature, T, except for ^(H^ |e _ ), which is in terms of the electron temperature, T e . 
However, it is expected that electrons and H 2 should be thermalized to the gas kinetic 
temperature via collisions, so we make no distinction between T e and T in our calculations. 
tThis expression includes contributions from the state-specific reactions H + + ( 3 Pj)0, 
where J = 0, 1, and 2, and makes assumptions regarding the relative populations in these 
fine-structure levels. Rate coefficients for the state-specific reactions were computed by 



Stancil et al.1 fjl999[K 

References : (1) iTheard Huntressl dl974h: ( 2 ) iMcCall et al.1 d2004h ; (3) IWoodall et al 



2007); (4) Stancil et al 



1999}; ffilSmith et al.1 (11978ft: (6) 



1990); (8) iRakshit fc Warneckl lll98(j ): (9) [Rosen et al.1 ((2000|); (10) I.Tensen et al.1 (j200Q) 



Jones et al 



(198l|klZ) 



Mitchell 
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Table 2:: IR Absorption Line Parameters 







f LSR 


FWHM 


W x 


a{W x ) 


N(J, K) 


a{N) 


Molecule 


Transition 


(km s -1 ) 


(km s _1 ) 


(10" 6 fim) 


(10" 6 urn) 


(10 13 cm" 2 ) 


(10 13 cm" 2 ) 


H 3 + 


R(l,l) u 


4.9 


11.5 


4.98 


0.64 


20.7 


2.65 


H 3 + 


R(1,0) 


5.4 


10.0 


5.01 


0.69 


12.7 


1.75 


H 3 + 


R(l,l) 1 


4.8 


10.1 


3.49 


0.14 


16.0 


0.64 


H 3 + 


R(l,l) u 


50.6 


12.8 


5.52 


1.76 


22.9 


7.30 


H 3 + 


R(1,0) 


48.5 


11.5 


6.51 


1.69 


16.5 


4.28 


H 3 + 


R(i,i) 1 


47.7 


13.9 


5.57 


0.17 


25.6 


0.78 


H 3 + 


R(l,l)» 


63.5 


10.0 


4.82 


1.33 


20.0 


5.52 


H 3 + 


R(i,o) 


63.2 


11.0 


5.93 


0.51 


15.0 


1.29 


H 3 + 


R(i,iy 


64.7 


6.7 


2.92 


0.12 


13.4 


0.55 






VLSR 


FWHM 


W x 


a(W x ) 


iV(J) t hin 


o-(-Nthin) 


Molecule 


Transition 


(km s -1 ) 


(km s _1 ) 


(10" 5 nm) 


(10" 5 fim) 


(10 15 cm" 2 ) 


(10 15 cm" 2 ) 


CO 


R(0) 


4.9 


10.5 


3.33 


0.17 


1.52 


0.08 


CO 


R(l) 


4.7 


8.7 


1.90 


0.28 


1.30 


0.19 


CO 


R(2) 








0.09 




0.07 


CO 


R(3) 








0.05 




0.04 



Notes: Columns 3 and 4 are the line-center velocity and velocity full-width at half- 
maximum (including instrumental broadening effects) found by a Gaussian fit to the ab- 
sorption feature. Columns 5 and 6 are the equivalent width, W\, and its la uncertainty. 
Columns 7 and 8 are the column density in the lower state, N, and its la uncertainty. The 
i?(l,0) line from the 5 km s _1 cloud and the R(l, l) u line from the 50 km s -1 cloud are 
severely blended. Note the difference in units for the equivalent widths and column densities 
of versus CO. In calculating a(W\) for the R{2) and R(3) lines a FWHM of 10 km s _1 
is adopted. Values reported for the R(2) and i?(3) lines of CO are la upper limits on the 
equivalent width and column density in the 5 km s _1 component. CO column densities are 
calculated assuming optically thin conditions, and are likely lower limits. For comparison, 
a curve-of-growth analysis with 6 = 2 km s _1 results in iV(0) = 1.97 x 10 15 cm~ 2 and 
N(l) = 1.49 x 10 15 cm" 2 . 
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Table 3:: Results from OH + and H 2 + Spectra 



^LSR 


iV(On + ) 


a[N(OR + )\ 


7\ T / T T / \ _|_ \ 

iV(o-H 2 + ) 


a[iV(o-H 2 + )] 




(km s _1 ) 


(10 12 cm" 2 ) 


(10 12 cm" 2 ) 


(10 12 cm~ 2 ) 


(10 12 cm- 2 ) 


/h 2 


0-11 


29.7 


1.30 


4.57 


0.72 


0.04 


11-17 


11.7 


0.55 


0.77 


0.31 


0.02 


17-21 


7.93 


0.36 


0.15 


0.15 


0.01 


21-33 


18.7 


1.00 




0.64 


< 0.03 


33-42 


8.94 


0.56 


1.99 


0.53 


0.08 


42-55 


54.4 


2.75 


7.26 


0.92 


0.04 


55-75 


34.7 


1.70 


7.12 


1.25 


0.07 



Notes: Velocity intervals were chosen to roughly correspond to the larger absorption 
components in the fit to the OH + and H 2 + spectra. /h 2 is calculated using equation 
(J5J) and assuming an ortho-to-para ratio of 3 for H 2 + . In the 21-33 km s -1 component 
the 3(j uncertainty in the H 2 + column density is used to determine the upper limit 
on /h 2 . Observations of ort/io-H 2 + toward W51 performed by the WISH (Water In 
Star-forming regions with Herschet) key program show absorption at 22.5 km s -1 at 
a position 60" away from the PRISMAS pointing (17" fro m the IR pointing). T he 



0.38 x 10 12 cm" 2 



flWvrowski et al.ll2010h 



inferred column density is A^(o-H 2 + 
consistent with our upper limit — and, when taken with our OH + column density, gives 
a m olecular hydrogen fract ion of 0.02 in that component. It should also be noted 
that IWyrowski et al.l (120101 ) find a column density of iV(o-H 2 + ) = 4.5 x 10 12 cm -2 
in a component at 6 km s _1 , in excellent agreement with our findings, validating 
our assumption that the slightly different sight lines targeted by the THz and IR 
observations should indeed probe similar material in the nearby diffuse cloud. 
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Table 4:: Molecular and Atomic Abundances in the Dif- 
fuse Cloud Toward W51 





N(X) 


W LSR 


FWHM 




Species 


(cm -2 ) 


(km s _1 ) 


(km s" 1 ) 


References 


H 


(1.39 ±0.3) x 10 21 


6.2 


5.6 


1 


H 2 


(1.06 ± 0.52) x 10 21 


3-10 




2 


H 3 + 


(2.89 ±0.37) x 10 14 


5.0 


10.5 


3 


OH+ 


(2.97 ±0.13) x 10 13 


0-11 




3 


H 2 + 


(6.09 ±0.96) x 10 12 


0-11 




3 


C+ 


(4.0 ±0.4) x 10 17 


0-11 




4 


CH 


(3.7 ±0.2) x 10 13 


3-10 




5 


CO 


(2.81 ±0.21) x 10 15 


4.8 


9.6 


3 



Notes: Values of i>lsr and FWHM describe the range of velocities over which the 
corresponding column densities were determined. In cases where vlsr and FWHM 
are each given by a single value, they represent the line center and full width at 
half maximum for a Gaussian fit to the absorption feature. In cases where vlsr is 
given by a range, the column density was determined by integrating between these 
two velocities. The H 2 + column density was determined from our observations of 
ort/io-H 2 + and an assumed ortho-topara ratio of 3. 

References : (1) Calcu lated from the H I absorption line parameters for G49.5-0.4e 
tabulated in IkooI ( 1997 ). but rescaled for a sp in temperature of 100 K; (2) estimated 



from iV(CH) reported by iGerin et al.l (12 01 Oaf) and the empiric al relat i onshi p between 



CH and H 2 repo rted bvf sheffer et all (EoS); (3) this work; (4) iGerinl J2012i in press); 
(5) lGerin etaD feplOah 



